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We argue for a "parametrized post-Friedmanian" approach to linear cosmology, where the history of 
expansion and perturbation growth is measured without assuming that the Einstein Field Equations 
hold. As an illustration, a model-independent analysis of 92 type la supernovae demonstrates that 
the curve giving the expansion history has the wrong shape to be explained without some form 
of dark energy or modified gravity. We discuss how upcoming lensing, galaxy clustering, cosmic 
microwave background and Lya forest observations can be combined to pursue this program, and 
forecast the accuracy that the proposed SNAP satellite can attain. 



I. INTRODUCTION 

Modern cosmology is in a somewhat equivocal state 
of affairs. The good news is that a recent avalanche of 
high-quality data are well fit by an emerging "standard 
model" whose roughly ten free parameters are being con- 
strained with increasing precision [1-11]. The bad news 
is that this emerging model is more complicated than an- 
ticipated. There is not one kind of invisible substance, 
but three: cold dark matter (CDM) to explain cluster- 
ing, dark energy to close the Universe and a smidgeon of 
massive neutrinos [12] that may well be too small to be 
cosmologically important. Moreover, problems involving 
small-scale clustering have triggered increasingly com- 
plicated models for the dark matter — self-interacting 
CDM [13-22], annihilating CDM [23], warm dark mat- 
ter [24-29,13,16,17], fuzzy dark matter [30] and fluid 
dark matter [31], to mention a few — and a large num- 
ber of dark energy models {e.g. [32,33]) have appeared 
where Einstein's single parameter A is replaced by a 
"quintessence" field that varies temporally and perhaps 
spatially. 

This perceived profusion of bells and whistles has 
caused unease among some cosmologists [34-38] and 
prompted concern that these complicated dark matter 
flavors constitute a modern form of epicycles. There have 
been numerous suggestions that, just as in the days of 
Ptolemy, the apparent complications can be eliminated 
by modifying the laws of gravity [39-49] For instance, 
there are scalar-tensor theories that can reproduce the 
observed accelerating Universe without invoking a cos- 
mological constant [43-45]. 

It will undoubtedly take time to settle these issues. 
In the interim, however, it is desirable to quote cosmo- 
logical measurements in a language that is fairly theory- 
independent. This is the topic of the present paper, de- 
limited to gravity in the linear regime, on large scales. 
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FIG. 1. Solid curve shows a currently popular model 
for the evolution of the effective cosmic mean density 
peff(-z) oc H(z) 2 . This curve uniquely characterizes the space- 
time metric to zeroth order. The horizontal bars indicate 
the rough redshift ranges over which the various cosmological 
probes discussed are expected to constrain this function. If 
the Friedmann equation is correct, then p e fi(z) = p(z). Since 
the redshift scalings of all density contributions except that 
of dark energy are believed to be straight lines with known 
slopes in this plot (power laws), combining into a a simple 
quartic polynomial, an estimate of the dark energy density 
px(z) can be readily extracted from this curve. Specifically, 
p oc (1 + z) 4 for the Cosmic Microwave Background (CMB), 
p tx (1 + z) 3 for nonrelativistic matter (like baryons and 
CDM), p oc (1 + z) 2 for spatial curvature, p oc (1 + z)° for a 
cosmological constant and p oc (1 + z ) 3 ^ 1+w ^ for quintessence 
with a constant equation of state w. Error bars are for our 
SNAP SN la simulation. 

A. Parametrizing our ignorance with two functions 

There is general consensus that space-time is well de- 
scribed by a metric, at least on non-microscopic scales, so 
competing theories of gravity differ in their predictions 
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for how the metric evolves with time and responds to the 
presence of matter. General Relativity (GR) makes spe- 
cific predictions in the form of the Einstein Field Equa- 
tions. The parametrized post-Newtonian approximation 
[50,51] has been highly successful in describing possi- 
ble departures from GR in the solar system, and the 
parametrized post-Keplerian approximation [50,51] has 
done the same for for binary orbits, notably binary pul- 
sars. Although GR has so far passed all experimental 
tests with flying colors, these tests tend to probe only 
the present epoch. Broad classes of theories have been 
shown to evolve into GR at late times, so substantial 
departures from the Einstein field equations are still pos- 
sible in cosmology. 

Assuming merely the cosmological principle (that our 
position is not special), it follows directly from the ob- 
served near-isotropy of the cosmic microwave background 
(CMB) that the Universe is approximately isotropic and 
homogeneous. This means that the large-scale metric 
consists of small perturbations on a Friedman- Robertson- 
Walker (FRW) background [52]: 
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for some function a(t) that describes the expansion his- 
tory of the Universe, k = —1, or 1 for an open, flat 
and closed Universe, respectively. If GR is correct, then 
a evolves according to the Friedman equation [52] 
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where we define the mean matter density p(z) to include 
a curvature contribution pu = — 8 ^ k Q a i ■ Other theories 
make different predictions. As summarized in Figure 1 
and reviewed in the next section, the function H(z) is 
directly measurable is a variety of ways that make no as- 
sumptions about dark matter properties or gravitational 
field equations, but rely on spacetime geometry alone. 

Let us now turn to the evolution of perturbations. The 
evolution equations are by definition linear to first or- 
der, regardless what underlying theory of gravity is be- 
ing linearized. Assuming that this theory is well approx- 
imated by partial differential equations on large scales, 
these decouple into ordinary differential equations for 
each Fourier mode, since we are perturbing around a ho- 
mogeneous metric and derivatives become local in Fourier 
space. In short, the key observable predicted by theory 
is the linear growth factor g(z,k), conventionally defined 
so that a density perturbation S(k) with wave vector k 
evolves with redshift as 



5( z ,k) = ^^«y(0;fc). 



(3) 



A wide variety of upcoming observations will probe 
g(z,k), including galaxy and Lyman a forest clustering, 



the cosmic microwave background, gravitational lensing 
and type la supernovae. When multiple matter com- 
ponents are present, there will of course be one growth 
factor for each one, whose evolution couple. For photons 
and neutrinos, the full phase-space distribution is rele- 
vant to the dynamics. On large scales, a broad class of 
models make the simple prediction that g(z, k) is inde- 
pendent of k. 



B. Is matter the matter? 

It is obviously premature to guess as to whether up- 
coming data will favor dark matter or modified gravity. 
However, it is worth noting that the distinction between 
these two cases is in fact somewhat blurry. As pointed 
out by Eddington [53], we can always choose to define 
matter as that which equals (8irG)~ 1 times the Ein- 
stein tensor, in which case the Einstein field equations 
G>„ = QnGT^v become satisfied by definition. This 
"matter" T^ v would have very strange properties for a 
generic metric, so the predictive power of GR arises from 
the fact that observed matter has simple, formalizablc 
properties. 

For many modifications of GR, the so defined T^ v is 
in fact rather regular. Indeed, in the familiar case of 
adding a cosmological constant, it is largely a matter 
of taste whether to call it modified gravity or dark en- 
ergy, corresponding merely to whether we insert it on 
the left or right hand side of the Einstein equations. In 
the currently popular class of GR generalizations known 
as scalar-tensor theories [52,50,44], of which Brans-Dicke 
theory is a special case, we can play the same trick: the 
Friedman equation (2) remains valid in the conventional 
(Jordan) frame [44] if we augment the density p with the 
extra "dark energy" term 
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where the theory is specified by the free functions F($), 
£/($) and Z($) of the dilaton field $, p m is the density 
of conventional matter, and 
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is the factor by which G e ff, the effective small-scale New- 
ton's constant measured in laboratory experiments today, 
differs from the "bare" one Gbarc which is truly constant. 

As another example, consider a flat toy model where 
equation (2) is replaced by 



H 2 = SnGp 2 
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for some constant p», i.e., a model where p cx H rather 
than p cx H 2 . An observer assuming that the Friedman 
equation equation (2) is valid would then conclude that 
there is a "dark energy" density p x given by 
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P=^-p=^(l + zf- Po (l + z) 3 1 (7) 
p* p* 



a quantity which goes negative for small p. Here po de- 
notes the current density of baryonic and dark matter. 

The modified behavior of linear perturbations in non- 
GR theories can clearly be attributed to new "matter" 
as well. Indeed, it has recently been shown [54] that the 
linearized properties of very general types of dark matter 
can be conveniently characterized by an effective sound 
speed and an effective viscosity. 

In light of this theory/matter ambiguity, we propose to 
use the Friedman equation as a definition of an effective 
mean density 
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Thus p e ff (z) = p(z), the true density, if the Friedmann 
equation is correct. As illustrated in Figure I, expressing 
observational results in terms of p e ff (z) rather than H(z) 
has the advantage of directly visualizing the dark matter 
problem. The contribution to p(z) from vacuum density 
is of course constant. Curvature density scales as (1 + z) 2 , 
matter density as (1 + z) 3 , radiation density as (1 + z) 4 , 
and a "quintessence" . component with constant equation 
of state w as (1 + z) 3(1+w \ 

A large body of work has been performed during the 
last few years on constraints and forecasts for more spe- 
cific models. Cast in our present notation, a particu- 
larly useful and popular parametrization of p(z) has been 
[55-74] 

p(z) = p° m (l + zf + p°(l + zf +p° x (l + z) 3(1+ ^ (9) 

for z <C 10 3 , where p {) m , p° and p\ are the present den- 
sities of matter, curvature and dark energy, respectively, 
and the dark energy equation of state w may in turn vary 
with rcdshift. In this paper, we treat p c s(z) rather than 
w(z) as the free function for two reasons: 

1 . It can be estimated from the data without assump- 
tions about gravity (the Einstein field equations) or 
about the dark energy (which may be out of equi- 
librium and not have any simple equation of state 
uniquely determining its pressure from its density). 

2. It is more directly related to observed data since, 
as will be elaborated below, it is given by the first 
rather than the second derivative of a measured 
function. 

A similar choice was made in [75] as well. 



Much progress has also been made on a reconstruc- 
tion program, aiming to reconstruct the microphysics of 
specific theories from data. The present work comple- 
ments this approach, since such reconstructions of, e.g., 
the quintessence potential [56,64,61,67,68] or the func- 
tions specifying scalar-tensor gravity [43] can be per- 
formed with Pcs(z) and g(z, k) as a starting point. 

The rest of this paper is focused on the Oth order func- 
tion p ff (z), with the 1st order function g(z,k) studied 
in a companion paper [76]. In the next section, we re- 
view techniques for measuring p c g(z), compute the con- 
straints from current SN la data and make forecasts for 
the proposed SNAP satellite. We discuss our conclusions 
in Section III. 



II. MEASURING p eH (z) 

The bars at the top of Figure I show the redshift ranges 
over which various types of observations are sensitive to 
Pes(z). We will first review them briefly in turn, then 
return the SN la in greater detail. 



A. Space-time geometric tests 

All cosmological tests based on luminosity, angular 
size and age (see e.g. [52,77]) can be described as noisy 
measurements of some quantities x n at redshifts z n , 
n = 1, ...,N, modeled as 



x n = a\nd(z) + b n + s n , 



(10) 



where a and b are constants independent of p e s(z), the 
function d incorporates the effects of cosmology and e n 
is a random term with zero mean ((e„) = 0) including all 
sources of measurement error. As briefly review below, 
these observables all probe different weighted averages of 
the quantity p e ff(z) -1 / 2 . 

For luminosity tests like SN la, x n is the observed mag- 
nitude of the n th object and d is the dimensionless lumi- 
nosity distance [52]: 

d lum = (1 + z)^l, V (z) = H f Jf-, (11) 
K Jo n \ z ) 

where the k = yjl — flo| is the spatial curvature, and 
51o = 87rGp o ff(0)/3-ff 2 is the fraction of critical density 
contributed by p c ff today. S(x) = sinhx, x or sinx when 
k = — 1, or 1, respectively. From the definition of mag- 
nitudes, a = 5/ In 10 and b n = 25 + 5 lg[if^ 1 c/lMpc] + 
M n , where M n is the (hopefully known) absolute mag- 
nitude of the n th object. The errors e n include errors in 
extinction correction and intrinsic scatter in the "stan- 
dard candle" luminosity. 

For tests involving the observed angular sizes 6 n of ob- 
jects at redshifts z\, zn with (hopefully known) linear 
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sizes Di, ...,Dn, we can take d to be the dimensionless 
angular size distance [52]: d ang — di um /(l + z) 2 . This 
gives x n = \n9 n (angles measured in radians), a == — 1 
and b n = ln[HoD n /c]. For such tests (e.g. [78,79]), e n 
includes scatter in the "standard yardstick" size. 

For tests involving estimates t n of the age of the Uni- 
verse at redshifts z n , we define x n = In Hot n . Setting 
a = 1 and b n = 0, this gives [52] 

f 00 dz' 

d ^ = H °J z { i + z >)H(z>y (12) 

Equation (10) implicitly assumes that the the mea- 
surement errors e n for additive. This is a popu- 
lar approximation in many cases, notably for luminosity 
distance where the actual measurement is a magnitude 
proportional to In d, and tends to be accurate if the rela- 
tive error (on flux, angular size, age, etc.) is small (<C 1). 
In other cases, a more appropriate error model should be 
employed. 

Ideally, we would like to measure the function p e e(z) 
in ways that are sensitive only to the global spacetime 
geometry and do not depend on g(z,k). The measure- 
ments that come closest to this ideal are arguably those 
involving the brightness, angular size and age of distant 
objects, in situations where the calibration of the "stan- 
dard candle", "standard ruler" or "standard clock" de- 
pends mainly on microphysics. However, great care must 
be taken when constraining competing theories of grav- 
ity etc. (as opposed to more mundane dark energy), since 
slight time variations in microphysical constants could 
affect the standard calibrators and masquerade as modi- 
fications to p e g (z). For instance, time- variation in a non- 
gravitational "constant" (say the fine-structure constant 
"alpha") could alter the intrinsic energetics and time- 
evolution of a supernova explosion, making it a bad stan- 
dard candle. If such an effect were present but neglected, 
then we would draw an incorrect conclusion about the lu- 
minosity distance to the supernova and consequently also 
about the cosmic geometry. 

Terrestrial constraints are are now so strong that time- 
variation of the fine structure constant a is unlikely to 
cause such misinterpretations: variations in a are limited 
to around a part per million over cosmological timescales 
[80] (c.f. [81]). A more serious concern in this context is 
time variations in the effective gravitational constant, ex- 
pected in many scalar-tensor theories according to equa- 
tion (5). Here the observational constraints are weaker, 
permitting variations of a few percent over cosmologi- 
cal time scales [82,83], which could among other things 
cause distant supernovae to be systematically brighter or 
fainter than nearby ones. 

As Figure I indicates, the function p c s(z) is probed by 
the growth of clustering as well, via the power spectrum 
growth in the Lya forest, weak lensing and galaxy surveys 
as well as via the abundance evolution of galaxy clusters 



and lensing systems. However, these probes all depend on 
the function g(z, k) as well, i.e., on the next order of per- 
turbation theory. The Alcock-Paczynski test elegantly 
eliminates the dependence on g(z, k) caused by number 
density evolution, but involves certain bias-related issues 
[84,85,61]. 

Although the CMB power spectrum Ci depends on the 
growth g, the effect is "vertical" rather than "horizontal" . 
In other words, as long as some detectable acoustic peak 
structure is visibly preserved, it should be possible to ex- 
tract the angular scale corresponding to the horizon size 
at recombination fairly independently of g. With addi- 
tional knowledge of the dark matter density, this provides 
a clean measurement of di um at z ~ l0-\ 

If GR is correct so that g and the vertical structure 
of the acoustic peak heights can be computed from first 
principles, then the CMB is a sensitive probe of the 
high-redshift asymptote in of the p e s curve in Figure 2, 
where the Universe becomes matter dominated. Writing 
p m (z) = p^(l + z) 3 , the MAP and Planck satellites are 
expected to measure p", to accuracies of 10% and 1%. re- 
spectively. However, in contrast to, e.g., SN la, the CMB 
is mainly sensitive to the physical matter density p m , not 
to the plotted quantity p m /h 2 , since this is what matters 
for the acoustic oscillations at z > 10 3 . The ultimate ac- 
curacy with which we can measure p^i/h 2 (equivalently 
Q m ) is therefore likely to be limited not by CMB issues 
but by our ability to measure h, either directly or by 
combining various complementary measurements [86]. 

B. p e ff(z) from existing SN la data 

Of the various low-redshift probes of p e e(z) that we 
discussed, type la supernovae are arguably the most sen- 
sitive at the present time. We will therefore compute 
constraints on p e s(z) from current SN la data and make 
forecasts for the accuracy attainable with future mea- 
surements. 

Figure 2 shows the measured luminosity distance for 
92 SN la, combining the high redshift sample of 42 
from [87] with the 50 reported in [88] with the MLCS 
method. Apart from systematic and calibration issues, 
our desired function p e s(z) is simply given by the deriva- 
tive of the curve around which they scatter: p c s(z) — 
3HQ/8wG e gri'(z) 2 oc h 2 /r]'(z) 2 .* Since the slope is seen 



*From here on, we will limit our discussion to the cur- 
rently favored case of flat space, i.e., k — 0, where r)(z) = 
dium{z)/(l + z) is directly measured for each supernova. 
In the general case, r](z) = k~ L sm~ 1 [kdi um (z) / (1 + z)] or 
r](z) — k^ 1 sinh -1 [kdi um (z)/(l + z)\ depending on the sign of 
the curvature, introducing an extra degeneracy with k that 
has been exhaustively treated elsewhere [89-91,86]. Since 
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to fall, the cosmic density clearly increases at high red- 
shift. If standard GR is correct (p e g = p), then this fact 
that dp e H jdz > corresponds to the weak energy condi- 
tion [92,75] that the density measured by any observer is 
non-negative, and is equivalent to the pressure-to-density 
ratio constraint p/p > — 1. A more detailed reconstruc- 
tion of p{z) is shown in Figure 3, and was computed by 
fitting a homogeneous quartic polynomial (dashed curve 
in Figure 2) to r)(z). Using the notation of equation (10), 
the fit was done by minimizing 



71=1 



x n - a\nd(z n ) - b n 



(13) 



where a n are the quoted 1 sigma magnitude errors 
[87,88]. This gives the fitting parameters (polynomial co- 
efficients) as linear combinations of the data x n . The 1 
sigma errors in Figure 3 are obtained by propagating the 
errors on these fitting parameters.^ The figure clearly 
illustrates the familiar fact that dark energy/modified 
gravity is required: the slope of the /9 e ff(z)-curve is too 
shallow at low redshift to be explainable in terms of ordi- 
nary matter (p oc [l + z] 3 ), spatial curvature (p oc [1 + z] 2 ) 
or some combination of the two. 

Since p c s(z) oc h 2 /rj'(z) 2 7 the SN la give us a direct 
measurement of p c if(z)/h 2 , where h is the Hubble param- 
eter today in units of 100 km/s/Mpc. Changing h thus 
shifts p c g vertically on a logarithmic plot like Figure 2 
without affecting its problematic slope. 

To assess the sensitivity of our results to method de- 
tails, we show alternative analyses in figures 4 and 5. 
Here we replaced the polynominal fitting function by 
piecewise quadratic and piecewise linear interpolating 
functions for rj(z), respectively. For the quadratic case, 
we enforced that the function have continuous derivative, 
thereby making p c ff a continuous function. The results 
are seen to be quite robust, in all cases agreeing well with 



CMB experiments are measuring the acoustic peak locations 
to increasing precision, and these are mainly sensitive to cur- 
vature (given by p c ff (0) /h 2 ) , degeneracies between different 
evolution scenarios for p e s(z) will be more important than 
the curvature degeneracy in the post-MAP era. 

^Equation (13) assumes that the magnitude errors on the 
different supernovae are uncorrelated. It is straightforward 
to generalize the procedure to the case of correlated errors, 
where the NxN magnitude covariance matrix is non-diagonal 
given a physical model for the nature of these correlations. 
A perfectly correlated term, corresponding to a common off- 
set in all magnitudes, would have no effect on the results 
since it is only the derivative that matters in Figure 2, but 
redshift-dependent offsets can of course bias the results. In- 
deed, a deeper understanding of possible systematic errors of 
this type is the main challenges for taking full advantage of 
the statistical power of future experiments like SNAP. 



the "concordance" model and indicating a shallow slope 
at low redshifts incompatible with just matter and cur- 
vature alone. A very nice non-parametric SN la analysis 
was recently performed by Wang & Garnavich [75] , and it 
is interesting to compare the two studies since they were 
performed concurrently and independently. The quan- 
titative agreement is good taking into account that [75] 
treat the dark energy density p{z) — /o^(l + z) 3 rather 
than p c fi{z) as the free function and marginalize over p m 
imposing the weak energy condition. The SN la data 
thus speak loud and clear: the dark energy puzzle is ro- 
bust to theoretical assumptions and method details. 
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FIG. 2. Crosses show the fuminosity distance for 92 SN 
la. From top to bottom, solid curves correspond to models 
(fi m ,f2 A ) = (0,1), (0.38,0.62) and (1,0), respectively. The 
middle curve is almost indistinguishable from the best fit 
quartic polynomial (dashed). The density history p c ff(z) is 
simply the squared inverse slope of this curve. Scatter in- 
creases with z since the relative distance errors are roughly 
constant. 
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FIG. 3. Zoom of Figure I showing constraints on p e s(z) 
from actual and simulated data, assuming a flat Universe. 
Solid black curve shows best fit to the 92 SNla, corresponding 
to the polynomial fit shown in Figure 2, and yellow/light grey 
area shows the associated 68% confidence region. Green/dark 
grey area shows the corresponding 68% confidence region from 
our SNAP simulation, for a fiducial model with f2 m = 0.38, 
Qa = 0.62 (red/grey curve) whose two components are shown 
as dashed lines. Error bars are for the non-parametric re- 
construction of Section II C, spaced so that measurements of 
neighboring bands are uncorrelated, and are identical to those 
shown in Figure I. 
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FIG. 4. Same as previous figure, but using a quadratic 
spline instead of a polynomial to fit r)(z). Bin widths were 
Az = 0.2 for the actual data and Az = 0.1 for SNAP. 
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FIG. 5. Same as previous figure but using linear interpola- 
tion to fit T)(z), giving a piecewise constant density. 



C. With upcoming SN la data 

The proposed SNAP satellite* would measure of order 
N = 2000 SN la out to redshifts around 1.2 or higher 
[61]. There has been substantial debate in the literature 
about the accuracy with which the equation of state can 
be reconstructed [60,61,64,66-75]. We will now compute 
how accurately SNAP can measure p e s(z). It is straight- 
forward to do this numerically in a "black box" fashion, 
and figures 3-5 show the result of doing this with sim- 
ulations in the same way as we did it for the currently 
existing data. However, because of the strong interest 
and substantial resources focused on designing such fu- 
ture probes, it is also worthwhile to gain intuition on 
how experimental specifications translate into model con- 
straints. We will therefore derive an analytic result which 
has the advantage of explicitly showing how the results 
scale when changing survey specifications, the smoothing 
scale, etc. 

1. An analytic error formula 

Defining f(z) = H(0)/H(z), we have rj(z) — 
Jo f{ z ')dz' and wish to know how accurately / can be 
recovered from noisy measurements of r\. In other words, 
we simply want to compute the derivative of a function 



'More information on the SNAP mission is available at 
http : / /snap. lbl.gov 
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measured with sparse sampling and noise. The errors on 
Pett(z) oc l/f(z) 2 then follow trivially from those on f(z). 

The accuracy with which f(z) can be measured is read- 
ily computed with the Fisher information matrix formal- 
ism [93,94] with infinitely many parameters to be esti- 
mated (the value of / at each z). The amount of infor- 
mation about these parameters in the SN la data is given 
by [90] 



where a = 5/ In 10, Am is the rms SN la magnitude er- 
ror, N is the number of supernovae, g(z) is their redshift 
distribution (normalized to integrate to unity), and 



w z '(z) 



dlnri(z) _ e(z-z') 



df(z') 



V(z) 



(15) 



since r)'(z) = f(z). Here 9 denotes the Heaviside step 
function (9(z) = 1 for z > 0, vanishing otherwise). Equa- 
tion (14) thus gives 

F(z, z') = J 9(z" - z)9(z" - z')h(z"Y x dz" (16) 

h(z"y x dz\ where (17) 



= / 

J m; 



' max(js,z') 

(Am) 2 ^) 2 
[ ' ~ a 2 N g{z) ■ 



(18) 



The covariance matrix C giving the best attainable error 
bars on our estimated parameters is the inverse of this 
infinite-dimensional matrix, C = F _1 . Fortunately, this 
inverse can be computed analytically, giving 

C(z, z') = -h'(z)5'(z - z') - h(z)S"(z - z'), (19) 

where 6 denotes the Dirac delta function. This follows 
from the identity 



/ 



C{z,z")F(z",z')dz" = 5{z-z') 



(20) 



which is proven by substituting equations (16) and (19), 
integrating by parts twice and using that 9'(x) — S(x). 
In practice we obviously want to smooth our estimated 
/-function (f(z), say) to reduce noise. Using weighted 
averages fi = J ifji(z)f(z)dz, the covariance matrix C%j 
for these averages is 

J C(z,z r )^(z)^(z r )dzdz' = J h(z)^' i (z)tp' j (z)dz, (21) 

employing equation (19) and integrating by parts. Once 
these errors on f(z) have been computed, it is of course 
trivial to map them into errors on p c e(z). When the 
relative errors are small (A/ <C /), we have simply 
Apeff/Poff ~ A///. 



2. Smoothing functions 
Using Gaussian smoothing functions 

1 -i(^r 



a 2 N r°° 

F(z',z") = 7 ^ g(z)w z ,(z)w z „(z)dz, (14) A/., 



55 J^Az 
the error bars Afi \ 
In 10 



e 2 ^ 



C jt ' are given by 



Am rj{zi) 



107rV4 (Az) 3 / 2 N 1 / 2 g(z i ) 1 / 2 



(22) 



(23) 



for the case where the smoothing scale Az is much 
smaller than the scale on which r\ and g vary substan- 
tially. In the same approximation, the dimensionless cor- 
relation coefficients between the different averages are 



(2 - n 2 )e ™ , where n = 



(24) 



is the number of smoothing lengths by which the two red- 
shifts are separated, so neighboring measurements are 
uncorrelated if they are separated by \/2Az. Another 
interesting case for our discussion below is that of "tri- 
angular" smoothing functions 



A(z) 



which gives 



ii\z- Zi \<Az, 
3 otherwise, 



Afi 



^lnlO 



Amrj(zi) 



(Az^N^giz,) 1 / 2 



(25) 



(26) 



in the same approximation. Here there are strong an- 
ticorrelations r = — 2~ x / 2 ss —0.7 between neighboring 
bins, and all longer-range correlations vanish. "Box- 
char" smoothing where we simply average in bins, tp(z) — 
9(Az/2 — \z — Zi\)/Az, may seem like a natural choice. 
This is a disaster, however, since the resulting 5-functions 
in V>i blow up when squared. This problem is easy to un- 
derstand physically. In a toy example where we have 
many SN la equispaced in redshift and estimate i]'(z) 
by simply differencing neighboring supernovae, averag- 
ing a segment of these derivative estimates will reduce to 
simply differencing the first and the last, placing all the 
statistical weight on merely two objects. 

As we will return to below, our black-box numerical 
method can also understood in terms of smoothing func- 
tions of a particular form. 

The exact result of equation (21) and the approxima- 
tion of equation (23) illustrate a number of issues rele- 
vant to measuring p e s(z) in practice, as discusses in the 
following subsections. 
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3. Redshift resolution 

Errors scale as (Az)~ 3 / 2 rather than (Az) -1 / 2 because 
of the derivative nature of what we are measuring, which 
means that very large number of SN la are needed to 
probe the small-scale structure of p e s(z). It also means 
that care must be taken in interpreting error forecasts to 
avoid coming away with a mislcadingly pessimistic im- 
pression. The superiority of SNAP over existing data 
is partly hidden in figures 4 and figures 5 since Az is 
halved for SNAP, roughly tripling error bars. The error 
bars from the numerical SNAP simulations in figures 4 
and figures 5 agree well with the analytic ones of equa- 
tion (25) with Az — 0.1 as well as those of [75], but 
are substantially larger than those from the polynomial 
fit shown in Figure 3. This is because the quartic poly- 
nomial enforces more smoothness, effectively increasing 
Az. The non-parametric error bars in Figure 3 from 
equation (23) using Gaussian smoothing are a factor of 
2 3 / 2 7r 1 / 4 rs, 4 smaller than those from equation (26) us- 
ing triangle smoothing, which is related to the fact that 
the triangle smoothing functions have standard devia- 
tions y/6 times narrower than the Gaussian ones. 

Which error bars are most appropriate depends on 
what the data are to be used for. Constraining mod- 
els where p e ff may change abruptly with z requires high 
redshift resolution and associated large error bars like in 
fugures 4 and 5. Most published models, however, pre- 
dict rather smooth functions p e s(z), so the small error 
bars from a low-order polynomial fit like in Figure 3 give 
the most accurate representation of how well SNAP could 
distinguish them from one another. 

4- High-z sensitivity 

Since errors scale as r)(z), one does worse at large z. 
However, since ij(z) flattens out around z ~ 1 and asymp- 
totes to a constant as z — > oo, this is not a fundamental 
limitation to probing very large redshifts — the challenge 
is simply to find suitable high-z standard candles or other 
reference objects to prevent g(z) from going to zero in the 
denominator. 

5. Optimal estimation 

Equation (21) also indicates what the optimal esti- 
mator is for non-parametric measurement of p e s(z). It 
is easy to verify that when the smoothing scale Az is 
much smaller than the scale h/\h'\ on which h varies sub- 
stantially, weighting the actual data with the function 
ip' i (z)h(z) and normalizing appropriately will recover the 
minimal Fisher matrix error bars. In practice, g(z) is of 
course not a smooth function but a sum of <5-functions 



at the observed SN la redshifts, so an additional correc- 
tion for the non-even spacing between adjacent SN la is 
necessary. 

The black box of % -minimization can also be under- 
stood in terms of smoothing functions and equation (21). 
Since x 2 depends quadratically on the free parameters, 
the best fit parameters and hence also the recovered func- 
tion f(z) depends linearly on the observed rj- values for 
each supernova, once the redshifts are given. In other 
words, we can write 

JV 

/i = 5>i(si)»fe (27) 

3=1 

for some function ipi , where rjj is the observed r\- value for 
the j th supernova. In an attempt to demystify the nu- 
merical results somewhat, we have plotted these weight 
functions in Figure 6 for estimating f(z) with SNAP 
at z ~ 0.8. We see that for both the quadratic spline and 
the linear interpolation cases, there is substantial ring- 
ing. Moreover, the measurement of 77' (0.8) is seen to in- 
volve supernovae over a very broad range of redshifts, not 
merely near z = 0.8. In contrast, the triangle smooth- 
ing function would be more local, giving ip'(z) = — 1 for 
z = 0.7 - 0.8, ip'(z) = 1 for z = 0.8 - 0.9 and zero 
elsewhere, simply estimating the derivative rj' (0.8) by 
subtracting measurements to the left from ones to the 
right. The Gaussian smoothing function behaves simi- 
larly, merely with less abrupt behavior. Just as sophis- 
ticated data analysis techniques developed by numer- 
ous groups have enhanced the science return of recent 
CMB missions, there may well be room for methodolog- 
ical improvement here. For instance, what is a good 
parametrization of / that leads to compact, easy-to- 
interpret weight functions? 
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ture of p e g(z). The reason that w(z) is so much harder to 
constrain accurately is that it involves a second derivative 
of the data, causing error bars to blow up as (Az)~ 5 / 2 . 



Redshift 



FIG. 6. Sample weight functions ip'(z) extracted from the 
X 2 -minimization are plotted for the quadratic spline (dashed) 
and linear interpolation (solid), both probing the derivative 
r]'(z) at z ~ 0.8. Excessive ringing drives up the error bars 
and the fact that the functions are so wide complicates the 
interpretation. 



III. DISCUSSION 

We have argued that the effective density evolution 
Pes(z) and the linear growth factor g(z, k) constitute the 
natural meeting point for theory to confront observation 
in terms of linear cosmology since they can be measured 
with very few assumptions about unseen matter and the 
underlying theory of gravity. This generalizes studies of 
the equation of state w. 

We analyzed current SN la data in this framework, 
and showed that the statistical constraints on p c g(z) 
are already fairly strong for z < 0.5. This alterna- 
tive approach confirms the results of previous SN la 
analyses [88,87,89,75], and provides a clean and model- 
independent case for either dark energy or modified grav- 
ity, since p e s(z) has the wrong logarithmic slope at low 
redshift to be explained by either matter (slope —3) or 
spatial curvature (slope —2). 

We estimated the accuracy with which the proposed 
SNAP satellite should be able to measure p D g(z), both 
with a numerical simulation involving polynomial fitting 
and by deriving a non-parametric analytic formula, and 
found that precision density measurements should be at- 
tainable in about a dozen independent redshift bins out 
to z ~ 1.2. Since errors grow as (Az)~ 3 ^ 2 as the bin 
width Az is reduced, the large numbers of SN la of SNAP 
are absolutely necessary to resolve the small-scale struc- 



Cosmic Microwave Backgroi 




Galaxy clustering 
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1 10 100 1000 

Comoving length scale A= 2u/k [Mpc] 

FIG. 7. Shaded regions show crude estimates of the ranges 
of scale and redshift over which various observations are likely 
to probe the growth factor g(k,z) over the next few years. 
The lower left region, delimited by the dashed line, is the 
non-linear regime where rms density fluctuations exceed unity 
for the "concordance" model from [9]. 

We focused on the 0th order function p e g(z), with the 
1st order function g(z,k) studied in detail in a compan- 
ion paper [76]. To place our discussion in its observa- 
tional context, Figure 7 shows the rough ranges of scale 
and redshift over which various observations are likely to 
probe g[k, z) over the next few years. Regardless of what 
detailed assumptions are made, it is clear from figures I 
and 7 that upcoming measurements of the CMB, the Lya 
Forest, weak lensing, cluster abundance and galaxy clus- 
tering have the potential to probe both functions back to 
high redshift over several orders of magnitude in scale, 
allowing numerous consistency checks between different 
types of measurements. Measuring g(z) for pressureless 
matter could provide an independent determination of 
H(z) [56]. 

It may be possible to extract further constraints on 
dark energy/modified gravity from going to higher or- 
der in perturbation theory or studying the nonlinear 
regime. Such nonlinear issues were indeed the motivation 
for some recent suggestions for dark matter and modified 
gravity [39,40,42]. Although many proposed quintessence 
fields do not cluster on small scales, early work in this di- 
rection [44,45] has discussed how temporal and perhaps 
spatial variations in Newton's constant may be an ob- 
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servable signature. 

In conclusion, the evidence for dark energy/modified 
gravity has emerged as one of the key puzzles of mod- 
ern cosmology. Fortunately there are at least half a 
dozen observational approaches that can directly probe 
the time-evolution of the effective mean density and clus- 
tering, so by going after this problem with the full arse- 
nal, there is real hope that the cosmology community can 
resolve it over the next decade. 

The author wishes to thank Arthur Kosowsky and Jim 
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Starobinski, Paul Stcinhardt, Mike Turner and Robert 
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